UR-1557 
ER/40685/926 
November 1998 

GLUON RADIATION IN TOP MASS RECONSTRUCTION: 
EFFECT OF HADRONIC W DECAYS 



L.H. ORR 

Department of Physics and Astronomy, University of Rochester, Rochester NY 14627-0171, USA 

E-mail: orr@pas.rochester.edu 



OO 

On 



T. ANDRE 

Department of Physics and Astronomy, University of Rochester, Rochester NY 14-627-0171, 
Department of Physics, University of Chicago, Chicago IL 60637, USA 
E-mail: troy@hep.uchicago.edu 



USA 



> 

o 
o 



> 
^1" 



T. STELZER 

Department of Physics, University of Illinois, Urbana, IL 61801, USA 
E-mail: stelzer@pobox. hep. uiuc. edu 

Top quark events in hadron collisions often contain additional hadronic jets from gluon bremsstrahlung off the quarks and gluons in 
the hard subprocesses. Such extra jets must be taken into account in attempts to reconstruct the momentum of the top quark from 
those of its decay products. We have performed a complete calculation of gluon radiation in top production and decay at the Fermilab 
Tevatron including hadronically decaying W bosons. In this talk we discuss the effect of gluon radiation on the reconstructed mass 
of the top quark, using various top mass reconstruction scenarios. Implications for the LHC are briefly discussed. rl 
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1 Introduction 

Measuring the top quark mass at hadron colliders re- 
quires reconstructing its momentum from its decay prod- 
ucts. Radiated gluons in top events can complicate the 
reconstruction process, because for example the jets from 
gluons can be indistinguishable from the jets in top de- 
cays. It is important to account correctly for these gluons 
because future top mass measurements will be dominated 
by systematic effects due to gluon radiation. 

Given a top event with an extra jet from a radiated 
gluon, what should we do with the extra jet? In partic- 
ular, should the extra jet be combined with the W and b 
quark to reconstruct m t ? The answer depends on where 
the gluon originated. If it was radiated from an initial 
state quark, then it is a correction to the production 
process that is not part of the top decay, and it should 
be ignored. If the gluon was radiated from one of the 
b quarks from the t or t decay, then it is itself part of 
the decay and should be included in the reconstruction. 
Suppose the gluon was radiated by the top quark itself - 
is it associated with top production or decay? In fact it 
can be either, depending on when the top quark went on 
shell. 

The point is that in a given event we cannot usually 
distinguish between the possibilities (even apart from the 
fact that they interfere), so we must consider top pro- 



duction and decay together in our treatment of gluon 
radiation. This has be£n done for_top production and 
decay at the Tevatron El and LHC □, without radiation 
from hadronic decays of the W bosons. But in the detec- 
tion modes in which at least one of the top quarks can be 
fully reconstructed from its decay products — the lepton 
+ jets and all-jets modes — one or both of the W bosons 
decays to quarks, which can themselves radiate. Radi- 
ation from hadronically decaying W's in top events was 
treated in the soft gluon approximation inu. The soft ap- 
proximation serves as a useful guide to the distribution 
of gluons and the relative importance of the various con- 
tributions, but it does not incorporate exact kinematics 
and cannot be used to study mass reconstruction. 

In this talk we present the results for an exact cal- 
culation of gluon radiation in top production and decay 
at the Tevatron with hadronic W decays fully taken into 
account. 



2 Gluon Distributions 

We have calculated the cross section for pp — > bbqqilvj 
and pp — > bbqq/qq/j where j is an extra radiated jet. 
This tree- level calculation is exact at 0(a^) and con- 
tains all spin correlations, top width effects, and inter- 
ferences. The center of mass energy is 1.8 TeV, and top 
and bottom masses are 175 and 5 GeV, respectively. He- 
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licity amplitudes are computed with the assistance of the 
MADGRAPH packagea. The results shown below are for 
the qq initial state that dominates in top production at 
the Tevatron; we have done the calculation for the gg 
and qg initial states El but do not show those results here. 

We apply the following kinematic cuts to all final- 
state jets (which in this parton-level calculation arc 
quarks and gluons) and to the charged lepton: 

Et 3 ,E T i > 15 GeV , 

\Vj\,h\ < 2-5 , 
AR jj: ARji > 0.4 . (1) 

These are meant to mimic experimental cuts, so that the 
partons are likely to appear in the detectors with enough 
angular separation to be distinguishable as separate par- 
ticles. They also protect the theoretical cross section 
from the soft and collinear singularities that appear at 
tree level. 

In the distributions we present below we will decom- 
pose the cross section into contributions from radiation 
associated with various parts of the process. These con- 
tributions are: 

i. Production-stage radiation, which comes from the ini- 
tial quarks or internal gluon line, or from the top (or 
antitop) quark before it goes on shell. 

ii. Decay-stage radiation that is part of the t or t decay; 
this is further subdivided into 

a. Decay-tb radiation from either of the &'s or from 
either of the t's after they go on shell. 

b. Decay- W radiation from the decay products of 
hadronically decaying W bosons. 

We make these distinctions in the parton-level calcu- 
lation based on kinematics to see how the various contri- 
butions behave; this cannot of course be done for a given 
event in the experiments. In principle the production- 
stage and decay-stage contributions can interfere with 
each other (with the exception of the decay-14 7 radia- 
tion, which cannot interfere with the other processes 
because the W is a color singlet). And although we 
do include all interferences in our calculation, in prac- 
tice the production-decay interference is very small for 
gluon energy thresholds large compared to the top width 
r t = 1.5 GeV, as in the present case. 

Figure ^ shows the distribution in pseudorapid- 
ity of the extra jet at the Tevatron for the lepton + 
jets case, i.e. for a single hadronically decaying W. 
The production-stage radiation, shown as a dashed his- 
togram, has the broadest distribution, populating most 
of the accessible rapidity range. The two decay-stage 
contributions are more centrally peaked. The decay-i& 
contribution (dotted histogram) is slightly larger, as it 
accounts for radiation from both the t and i, but the 
decay- W contribution (dot-dashed histogram) from a sin- 
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Figure 1: The extra jet pseudorapidity (jy) distribution (solid his- 
togram) and its decomposition in terms of production (dashed his- 
togram), decay-tb (dotted histogram) and decay- W (dot-dashed 
histogram) emission contributions, for the lepton + jets mode. 

gle W is similar in size and shape ±o the decay-i&, as was 
found in the soft approximation 13. The central region 
of the detector is populated by all three contributions, 
which means that distinguishing them will be challeng- 
ing at best. 

Results for the all-hadronic mode, where both W 
bosons decay to quarks, are similar, the main difference 
being that the decay- W contribution approximately dou- 
bles in size. 

The transverse energy distributions of the radiated 
jet are shown in Figure || for the lepton + jets (top) and 
all-hadronic modes. The spectra look quite similar for 
the various contributions, and extend to large values of 
transverse energy. 

3 Top Mass Reconstruction 

Given the difficulty of distinguishing extra jets radiated 
in the production stage from those from decay, we can 
ask what effect the extra jet has on mass reconstruction. 
First, however, we have to ask which jet is the extra one. 
It seems reasonable to assume that the gluon jet is the 
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Figure 2: The extra jet transverse energy distribution (solid his- 
togram) and its decomposition in terms of production (dashed his- 
togram), decay-tb (dotted histogram) and decay-W (dot-dashed 
histogram) emission contributions, for the lepton + jets mode (top 
plot) and all-hadronic (bottom plot) mode. 



one with the lowest Et, given the infrared singularity 
that characterizes emitted radiation. While it is true 
that of the final state jets, the gluon has the softest Et 
spectrum, the gluon has the lowest Et only just over 
half the time (just under half for the all-jets mode). Still, 
there is no method that is obviously better for identifying 
the extra jet. 

We reconstruct the top mass from the final state par- 
tons without assuming we know which jet is which, but 
we omit the jet with the lowest Et- For the lepton + 
jets mode we do the following. 
(1.) Drop the lowest Et jet. This leaves four jets. 
(2.) Find the jet pair with invariant mass closest to mw- 
(3.) Solve for the neutrino four-momentum using the 
charged lepton momentum and W mass constraint. 
(4.) Combine each of the W's with each of the remain- 
ing jets to give m t and mj. Choose the combination that 
minimizes the t-t mass difference. 

The mt distributions obtained from this procedure 
are shown in Figure ||. The top plot shows the distribu- 
tion for the top with the hadronically decaying W, and 
the bottom corresponds to leptonic W decay. In both 
cases we see a peak at the correct central value, where 
the procedure resulted in the correct mass. We also see 
smooth, reasonably flat high and low tails from wrong 
combinations, in addition to bumps in the low tails cor- 
responding to the omission of jets that were part of the 
decays. Note that these bumps appear in the contri- 
butions from decay-stage radiation, when all of the jets 
should be included in the mass reconstruction. Finally, 
the leptonically decaying W gives a sharper top mass dis- 
tribution because with no radiation from the W decay, 
there are fewer wrong combinations. 

The results shown in Figure [| are meant to be illus- 
trative and should not be taken as a direct representation 
of distributions measured in experiment. In particular, 
this calculation is at the parton level; we have not in- 
cluded backgrounds; and these distributions only include 
events with a radiated gluon. The effects of hadroniza- 
tion, energy resolution and detector effects, and back- 
ground will certainly make things worse. However there 
are certainly ways to improve on this simplistic analysis 
as well. For example, although we minimized mass differ- 
ences we did not cut on them explicitly. Figure |] shows 
the magnitude of the t-t mass difference on an event-by- 
event basis for the case shown above. It suggests that 
an absolute cut on the mass difference could reduce the 
tails in the mass distributions. Interestingly, ^-tagging, 
i.e. assuming we can identify b jets, does not improve 
the mass distributions much. This is because we do not 
include backgrounds. 

The results for the all-jets mode are similar, except 
that with both VF's decaying to quarks there are two 
more jets in each event, leading to more possible wrong 
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Figure 3: The reconstructed top mass distributions (solid histograms) in the lepton + jets mode for the hadronically decaying W (top plot) 
and the leptonically decaying W, and their decomposition in terms of production, decay-fft and decay-W emission contributions. 
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combinations and a corresponding increase in the tails. 
4 Conclusions 

In top quark physics, as statistics improve, systematic ef- 
fects associated with gluon radiation will dominate mea- 
surements of the top mass. We have added hadronic W 
decays to analyses of gluon radiation in top production 
and decay and presented some initial results here. We 
find that the contribution from radiation from a single 
hadronically decaying W is nearly as large as and com- 
parable in shape to the remaining decay-stage radiation 
from both the t and i. The presence of radiation from 
both the top production and decay stages complicates 
the reconstruction of the top momentum from its decay 
products and hence complicates the measurement of the 
top mass._Further analysis is in progress for the Tevatron 
and LHcli 
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Figure 4: The distribution in magnitude of the t-t mass difference 
for the events in the previous figure. 
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